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Abstract 

This paper considers Hamiltonian identification for a controllable quantum system with non-degenerate transi- 
tions and a known initial state. We assume to have at our disposal a single scalar control input and the population 

i 1 measure of only one state at an (arbitrarily large) final time T. We prove that the quantum dipole moment matrix 

i-C is locally observable in the following sense: for any two close but distinct dipole moment matrices, we construct 

discriminating controls giving two different measurements. Such discriminating controls are constructed to have 
three well defined temporal components, as inspired by Ramsey interferometry. This result suggests that what may 
appear at first to be very restrictive measurements are actually rich for identification, when combined with well 
designed discriminating controls, to uniquely identify the complete dipole moment of such systems. The assessment 
supports the employment of quantum control as a promising means to achieve high quality identification of a 
Hamiltonian. 

> 

I. Introduction 

fx. Quantum control has been receiving increasing attention [lj and one of its promising applications is 
to Hamiltonian identification pj by using the ability to actively control a quantum system as a means 
to gain information about the underlying Hamiltonian governing its dynamics. The underlying premise 

^ is that controls may be found which make the measurements not only robust to noise but also highly 

y—i sensitive to the unknown parameters in the Hamiltonian. Hence, although the performance of laboratory 
j> measurements may be constrained, the ability to control a quantum system has the prospect of turning 

JJp this data into a rich source of information on the system's Hamiltonian. 

In this paper, we consider the problem of identifying the dipole moment matrix (which is assumed to be 
^ real) of an N— level quantum system, initialized to a known state (ground state), from a single population 
measurement at one arbitrarily large time T. We suppose an ability to freely control the system with a 
time dependent electric field e(t). The measurements are obtained by (i) initializing the system's state at 
a known state |z), (ii) controlling the system with an electric field €k(t) and (iii) measuring the population 
of one state |/). This may be repeated for many controls (tk)k- A typical setting is for the system under 
study to be composed of a low pressure cloud of atoms or molecules with the interaction between the 
system and the control taking place on a time scale of ~ 100/s with the measurement of a population 
performed on a time scale of ~ 100ns. The large difference of time scales (lfs = I0~ 6 ns) motivates the 
assumption that the measurement is performed at an arbitrarily large time T. We prove that the dipole 
moment matrix is locally observable, that is, that for any two close but different dipoles, we can find a 
control which gives two distinct measurements to first order. 
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Identifying the dipole moment of a quantum system is important for fundamental problems in atomic 
physics such as the analysis of parity violation measurements [3|. It can also serve as a test bed for 
theoretical calculations using novel techniques such as relativistic many body perturbation theory Q. 
Various methods have been considered for measuring dipole moment matrix elements (see [|5j and 
references therein). The most common method is to measure the lifetime of a quantum state ||6), which 
is related to the norm of the dipole moment matrix element by Fermi's golden rule. Although this 
procedure has been extensively used to extract matrix elements between the ground and first excited 
states, it can be difficult to generalize the method to measuring matrix elements between two excited 
states. Problems such as branching and super-radiance Q complicate the interpretation of the data. Also, 
lifetime measurements only provide the norm of the matrix elements and give no information about the 
phases. Another commonly used method is multi-dimensional spectroscopy [[8j, which can reveal valuable 
information on the dipole coupling structure of the system (i.e, revealing which dipole elements are zero 
and non-zero). However, multi-dimensional spectroscopy is experimentally challenging to realize and 
often does not give quantitative information about the dipole moments. Other approaches have also been 
investigated including the measurement of the degree of induced polarization as a function of the detuning 
in a two photon process pj, stark shift measurements [9] and Rayleigh, Raman scattering (TO). 

In this paper, we consider the identification of the dipole moment matrix as an optimal control problem 
[11 1. We prove the existence of controls which make the identification from one population measurement 
a well posed problem (theorem [TJ). The controls we construct for this purpose are inspired from the two 
pulses separated by a large time interval used in Ramsey interferometry (see, e.g., JT2l page 148]). In 
Figure [TJ we illustrate this analogy and explain how to achieve, in three steps, a control to discriminate 
the real dipole moment between states \l) and \k). The field structure has the following components: 

1) Steer the system from its known initial state to the state \l) with a pulse e%. 

2) Perform a large number of Rabi flops between states \l) and \k) with a resonant control e cos(uik(r — 
Ti)), n>T< r 2 . 

3) Achieve a control pulse e 2 inducing an overlap with the measured population; such overlap will 
depend on the Rabi angle produced during the second step; thus the measured population will be 
sensitive to the dipole value between \l) and \k). 

Such discriminating controls are exploited in proving the key lemma [TJ later in the paper. 

The perspective above combined with control theory is motivated by three practical arguments. First, 
measuring a state population at one time T is a technique which can have a very high signal to noise 
ratio (~ 100). Second, technological progress with spatial light modulators (SLM) permits generating a 
broad variety of controls in the laboratory. Third, spectral phase interferometry for direct electric-field 
reconstruction of ultrashort optical pulses (SPIDER) [13| has proven to be a reliable technique to measure 
ultra-short pulsed fields. Hence, we are able to design a variety of controls and measure the actual controls 
created in the laboratory, along with the output signal. This scenario is suitable for parameter estimation, 
which is a well known problem in nonlinear system identification depending strongly on a careful choice 
of discriminating control inputs. 

Although the high potential of employing quantum control for Hamiltonian identification is clear, a 
full theoretical foundation is yet to be established. This effort is referred to as bilinear due to the form 
of the Schrodinger equation with the dipole and wave function entering multiplicatively. A vast literature 
exists for parameter estimation of bilinear systems, but to the best of our knowledge, identifiability results 
mainly rely on continuous time measurements linear with respect to the state. In a large set of quantum 
systems only the expectation value of a Hermitian operator can be measured, which is highly nonlinear 
with respect to the dipole. Le Bris and al [ 14 1 prove the observability of the dipole moment matrix when 
the population of all states are measured over an arbitrarily large interval of time. Methods to reconstruct 
the dipole from the measured data were proposed using nonlinear observers p"5) , p6| . A different setting 
is considered in [17|, [18| where it is supposed that one can prepare and measure the system in a set 
of orthogonal states at various times, and the available data is the probability to measure the system in 
a certain state when it was prepared in another; Bayesian estimation is used to reconstruct the energy 
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levels, the damping constants and the dipole moment from the measured data. We consider here the less 
demanding case where the only available measurement is the population of one state at one arbitrarily 
large time, and the initial state is known and coincides with the ground state. 

The paper is organized as follows. In section [II] we state the main result in Theorem [Tj and section III 



gives the proof of the Theorem. The proof is based on a lemma which we prove in section IV Finally 
concluding remarks are presented in Section |V| 



II. Observability of the quantum dipole moment 

A. Problem setting 

We consider a quantum system in a pure state described by the wave function G S. Here S is the 
set of N dimensional complex vectors of unit norm. The system interacts with an electric field (the real 
control input) e G St for some T > with 



& r = C 1 ([0,7];R) = if: [0,T] R £ \f{x)\dx < oo 



(1) 



For a given control e we measure the population of the state |/) at time T denoted as Pi/(e). We denote 
by H the free Hamiltonian (Hermitian matrix) and by p the dipole moment operator, also a Hermitian 
matrix. The initial state is an eigenvector of H denoted We consider a semi-classical model for the 
light-matter interaction [19|, and the dynamics of are given by the Schrodinger equation: 

° ' ' ' (H - e{t)p) \m) 



W)) 
W)) 



I*) 



(2) 



For all T > 0, we suppose that we can create any field e G St and that we can measure Pif(e). For M 
fields {ei, .., €m} we can collect the measurements {Pj/(ei), .., i-»/(eAf)}. Through ([2]) is a function of 
[i and a functional of e, and when necessary this explicit dependence will be written as Pj/(e, p). The aim 
of this paper is to explore the feasibility of estimating the dipole moment matrix p from the measured data 
{Pif(ei), .., Pif(eM)} using well chosen controls {ei,..,ejvf}- Below, P i f(e,fx) refers to the measurement 
achieved on the real system using a control e, and for any p ^ /i, Pj/(e,/t) is the estimated measurement 
which is obtained by simulating system ([2]) with the control e and coupling /}. 



B. Main result 

For all k < N we denote \k) as the eigenvector of H with associated eigenvalue E k . Throughout the 
paper, all matrices are written in the basis (|1) , .., \N)). For all k, I < N we specify af k = \l) (k\ + \k) 
We define 

M = Spm{a l x k \k, I < n such that Tr (/kx^) ^0}, M — dim(M) (3) 

with Tr being the trace operation. When all non diagonal elements of fx are non-null, M = M^zil jf^ 
main result is the following: 

Theorem 1. Consider a real symmetric matrix fi with zero diagonal entries and a real diagonal matrix 
H with non-degenerate transitions. Suppose that the system state in §2§ is controllable. Then for any 
positive constant a there exists a time T > and M fields (ei, .., €m) G St M such that the cost function 

M 

J-.fieM^Yl ( p iffa>fi - M e k,ii)) 2 (4) 
k=i 

is locally a-convex around p.. 
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In appendix [A] we provide the definitions of controllability, a— convexity and a matrix with non- 
degenerate transitions. 

A direct consequence of Theorem [T] is the local observability of the dipole moment matrix: 

Corollary 1. Under the assumptions of Theorem [7J the dipole moment matrix is locally observable. 

Here and throughout the paper, the norm of matrix p, noted refers to the max norm. 

Proof: Take a > 0. Theorem [T] implies that there exists a time T > and M fields (ex,. .,6m) 
G £t M such that the cost function J is locally a-convex around p. Hence there 3r > such that for 
all p G M with \\p — p\\ < r and p ^ p, J(p) > 0, and hence there exists e G {ei, ..,e M } such that 
P if (e,p)-P if (e,p) ^ 0. 

■ 

Remark 1. 77ze /oca/ a— convexity is a property stronger than the mere possibility to identify the dipole 
matrix. It states that the distinction between a dipole candidate p and the true dipole p can be observed 
(through the measurements aggregated in J) to first order in the distance \\p — p\\. This first order 
dependence of the measurement Pi f with respect to the dipole p is addressed in more detail in lemma [7J 

Remark 2. The cost function defined in Theorem [7J is locally a convex on M.. This suggests that all 
non-null matrix elements of p can be reliably estimated with well chosen controls, but we cannot say the 
same for null matrix elements. In practice, when dealing with an atomic or molecular quantum system, it 
is generally well known which matrix elements are null by applying selection rules (see Ax in, I.e. 7 in 
/[20]/j. Also, other methods, such as multi-dimensional spectroscopy, may be used to reveal the non-zero 
structure of the dipole moment matrix. 



C. Robustness 

This paper will not analyze the influence of noise effects on the estimation of the dipole moment. Some 
issues relevant to noise have been treated elsewhere J2TJ. We recall the three main sources of uncertainties: 

1) Uncertainty in the control e(t): a predetermined control cannot be perfectly created in the laboratory, 
and the actually created control cannot be measured with arbitrary precision. 

2) Uncertainty in the measured signal P%f(e). 

3) Uncertainty in the model: the bilinear Schrodinger equation ([2]) is a simplified model not taking 
into account dissipation, collisions and higher order matter-field interactions. 

We consider a simple noise model in order to give an estimate of the robustness of employing quantum 
control for Hamiltonian identification. Suppose that instead of measuring P i f(ek,p), one actually measures 
Pif(ek, p)+Vk, where i/& is the realization of a stochastic variable v of (finite) variance variy). We can then 
build a cost function J no isy with these measurement outcomes. Then it can be proven that J no i sy will reach 

its minimum at a point p no i sy , which belongs to a ball centered around p and of radius R no ise ~ \J var ^ 
(with a defined in Theorem [T]). This suggests that if the variance of the noise is a bounded function of 
the measurement time T and the controls in St, then the estimation of the dipole moment matrix can be 
arbitrarily precise by picking controls such that a is arbitrarily large. 



III. Proof of the TheoremQ] 

A. Existence of discriminating controls 

Theorem [T] is a local result and in order to prove it we need to study the first order variation of the 
measurement Pif with respect to the dipole moment p for a given time T and control e. We denote 
p' = pjpti the normalized dimensionless dipole moment operator. We denote p' lk = Tr (p'cr 1 ^) and 
the partial derivative of Pj/(e) with respect to p' lk . 

To prove Theorem [T| consider a class of controls which can be configured in terms of the physical 
framework sketched in Figure [T] with a suggestive analogy from Ramsey interferometry. In experiments 
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where shaped pulses are used as controls generated with spatial light modulators J22| , we can further 
expect that a good pulse to identify a certain dipole moment ji lk between states \l) and \k) can be achieved 
using the three steps set out in Section |IJ In particular, theorem [T] is based on the following lemma: 

Lemma 1. Suppose that fi is real, symmetric and has only zeros on its diagonal and Hq is real, diagonal, 
with non-degenerate transitions. Suppose system ([2]) is controllable. Then for all (I, k) such that pik ^ 
there exists £o > such that for all £ e]0, £o[ there exists T > and e G St such that 

1} wr k ^ = i + o(i) 

2) V{m,n}^{l,k} °*L( e ) = 0{1). 

Where 0(1) is a term bounded with respect to £ g]0,£ [- 

Remark 3. Lemma [7J claims that we can always find a time T and a control e which makes the data 
Pi/(e) arbitrarily sensitive to one element of [i while its sensitivity to the other matrix elements stays 
bounded. This is an existence result. In the laboratory however we are constrained by the measurement 
time T which has to be much smaller than the excited states lifetimes and the collision time as well as 
by the pulse energy and bandwidth which can not be arbitrarily large. These circumstances may result in 
the actual sensitivity being less than the theoretical estimate. Besides, the model (|2]) itself would possibly 



have to be modified for very strong fields [23] introducing extra unknown parameters in the dynamics. 



B. Proof of Theorem [7J 

Proof: To each pair of integers (l p ,k p ), l p < k p such that Tr (j//crx kp j ^ we associate a unique 

index p G {1, M}, and we define o v x = <ri pfep along with pi = Tr (//V§). 

According to lemma [lj 3£ > such that VO < £ < £ 3 Ti,..,T M and (ei,..,eu) G &r x x •• x £t m 
such that: 

1) for all p G [1 : M] (e p ) = £ + 0(1) 

2 ) WpgH(g = 0(l) 

We take T = max(T 1; .., T M ) and for all k G {1, .., M} we extend the definition of from [0, T k ] to 
[0, T] by taking e h (t) = for all t e}T k ,T]. 
We will use J : M — > R defined by: 

M 
k=l 

and then find 



so that for all p G {1, M} : = ^ + 0(|) and when p ± p> g^(p) = 0(\). We have: 

VV(/i) = ^(/ + O(0) 

where V 2 J(/x) is the Hessian of J at fi and / is the identity matrix. The smallest eigenvalue of V 2 J(p) 
scales as ^(1 + 0(£)), hence by taking £ small enough it can be made larger than a thereby reaching 
the conclusion above. 
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IV. Proof of lemma [T] 

A. Sketch of the proof: 

We introduce the dimensionless time scale r = |||i7 ||t and also T = |||i/o||^. 

In order to prove lemma jlj we need to find a control e G £ which makes |^ (e) arbitrarily large while 

c. P ^ (e) stays bounded. The concept is based on splitting the time interval [0,T] into three pieces [0, Ti] 
[Ti, t 2 \ and [r 2 , TJ. Over [r 1; r 2 J, playing the role of the time interval between the two Ramsey pulses, we 
will use a field resonant with the / — k transition producing extensive Rabi flopping between states \l) 
and \k). However, in order for this Rabi flopping to contribute to increasing |^(e), we need to drive the 
known initial state \i) to a particular state at r = T\ (i.e, an analog of the first Ramsey pulse) and also 
drive the final state |/) to a particular state \ip 2 ) at r = r 2 (i.e, an analog of the second Ramsey pulse). 
By^ increasing r 2 — T\ (i.e an analog of the time interval between the two Ramsey pulses), we will see that 
-9-7^ (e) can be made arbitrarily large. Besides, H having non-degenerate transitions, all other transitions 

dP 

will be out of resonance and hence we will show that as r 2 — t\ gets larger, V{m,n} 7^ {/, k} a-r*-(e) 
stay bounded. The particular controls ei and e 2 over the time intervals [0, T\\ and [r 2 , TJ are used to drive 
the system to the respective desired states and |^ 2 ) (see Figure [T]). Since we suppose the system to 
be controllable over the compact S, we can find such controls which are defined over the bounded time 
intervals [0, Ti] and [t 2 ,T]. Notice that such discriminating fields e (plotted in Figure [T]) may be difficult 
to create in the laboratory due to their very specific spectrum (a peak at coik and a broad distribution at 
other frequencies). This form of control is used here to construct an observability proof only. We are 
currently investigating the means to generate discriminating fields which may be readily created in the 
laboratory. 



Before starting the proof we introduce some notation. For two times r, t' G [0, T] we define the 
propagator U(t',t) such that \ip{r')) = U(t',t) \ip(r)). Rewriting ([2]) for U(r, 0) we obtain: 

i£pU(r,0) = ^-r(H -e(r)ii)U(T,0) 

U(0,0) = I (5) 
Pif(e) = |(/|f/(T,0)|.)| 2 

B. Proof of lemma [7] 

Proof: The proof of lemma [T] has two parts I and II separately treated below. 
i) Part I: Take two times Ti, r 2 , < T\ < r 2 < T. We can write (for any complex z we denote by 
z its complex conjugate): 

Pif{t) = I (/I U(T, r 2 )U(r 2 , npin, 0) \i) | 2 = zz where z = (f\ U(T, t 2 )U(t 2 , nMn, 0) |i> . (6) 
The differential of Pij is SPif = 25R (5zz) where 3?() denotes the real part. 



5z = (/|5C/(T,t 2 )L'"(t 2 ,ti)L'"(ti,0) \i) 
+ (f\ U (T, t 2 )5U (r 2 , Ti)U (ri, 0) |z) 

+ (f\U(T,T 2 )U(T*,T 1 )8U{n,G)\i) 



We find (see appendix B-A) 

d 



(7) 



Ik 



1 1 1 1 r 
Ufa, n) = ^YTTuUfa, n) / e(r)C/t(r, n)a l x k U(r, n )dr 

\\ H 0\\ Jrx 



(8) 
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Fig. 1. A good control e has three components (inspired from Ramsey interferometry) to enable the identification of /iik- The field ei is 
defined over [0, n] (i.e, an analog of the first Ramsey pulse) to steer the known initial state \i) to \ipi) = \l): \l) = U(n,0) |z). The field £2 
is defined over [t2,T] (i.e, an analog of the second Ramsey pulse) is such that |/} = J7(T,ts) \%j)2) where \ip2) is collinear to ^ • 

The relative phase is such that |^>a) = J7(r2,Ti) f^^^) where the propagator (/(r2,Ti) corresponds, for a long interval T2 — Ti, to a 
large number of Rabi oscillations with the control ecos(oj i ' fe (r — n)) resonant with the |Z) o transition. 



Denote for any m, n = 1, M : = E ^ H ^ n an d consider the control defined on [t\, t 2 \. 

e (r)=ecos(u;; fc (r-r 1 )) (9) 

where e is a small strictly positive real parameter. Take £ = p^|- The only remaining degree of freedom 
in the control over [7i,r 2 ] is £, which can be taken arbitrarily small. We define H' = j^H and 
u'mn = ( m \ H' \m) — (n\ H' Q \n). Note that u' mn = —oj' nm . We now rewrite ^ and ([8]) for the control given 



in (|9]) on the time interval [r\ , r 2 ] : 

d 



tfrUfan) = (^-ecos^r-rO^O^rO (10) 
° U(r 2 ,T 1 )= i ^U(r 2 ,T 1 ) r cos(u' lk (r - r^^ny^U^T^dr (11) 



n 



The goal is to show that -rA-Ufa, T\) can be made arbitrarily "large" while -5-7— U(r 2 , r x ) stays bounded. 
Note that all the terms in the integrand of ( [11] ) are bounded, and a rough estimate of the norm of 

-A-U{t 2) ti) gives a quantity proportional to (r 2 — Ti)£. Hence, we take r 2 — ti = implying the need 

1 i — i 

to have expressions for U(r, r x ) over a time scale on the order of To this end we state lemma 2 which 

gives such an approximation. 
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Lemma 2. Consider equation (10). There exists a Hermitian matrix K and £o > such that for any 
£ G]0, £o[ an d any r G [t\,T\ + ^] we have: 

U(t, n) = e - tH ^-^e<^ k+ ^ K ^-^ + 0(0 , 
where 0(£) ts sac/z ?/za? |0(£) w bounded and does not depend on r. 

We continue with the proof of Lemma [TJ and will come back to Lemma [2] in Section IV-C 



Using the expression of U(t,ti) given in lemma [2j the integrand in ( [TT] ) is: 

cos(a;; fc (r-r 1 ))f/ t (r,r 1 )ai fe f/(r,r 1 ) = 

e -^4^-K^)( T -n) (cos^r - r 1 ))e^o(-n)^ e -^ ) (r-n) N ) ^^-k^Xt-t + O (0 



In order to compute (11), we need the following result: 



cos(u4(r -T X ))e 



1 Jft 

2^ 



(12) 



+ s cos ( 2t 4( r " n))** + TT sin(2w; fc (r - n ))a, 



where we denote = — z |/) (A;| + z (Z|. In ( fT2[ ), the terms oscillating at frequency 2u;[ fe independent 
of £ will only contribute to 0(£) in ([TT]). We now focus on the contribution o f the term with <j l £ in ( [TT] ) 

calling for e- l ^^ k+eK ^ T ~ T ^al k e<^ a - +eK ^ T - T1 \ We find: (see appendix B-B) Vr 



e -<^^ fc +€ 2 ^)(r-r 1 ) a ;fc el (^^ fc +? 2 ^)(T-T 1 ) 



Ik 



o(C) 



Introducing ( 13 ) into (11), we find: 



(13) 



d 



d ^\k 



From now on, we take r 2 = T\ + ^ and obtain: 



^rt/(r 2) n) = iU(t 2 , n)(^a[ k + 0(1)). 



(14) 



(15) 



We define \ipi) = \l) and \ip 2 ) = T3^( r 2) r i)(IO Since the system is controllable there exists a 

time 7"i and a field e\ G £ ri such that U(ti, 0) |z) = \ipx), and there exists a time T and a field e 2 defined 
over [r 2 ,T] such that W(T,t 2 ) \f) = \ip 2 ). 

Since the state space is compact (here it is a sphere), we know that if the system is controllable, it 
is controllable in bounded time, and with bounded controls (see Theorem 6.5 in p4j). We can therefore 
define a bounded set of times \T(ipi, ipf)\(ipi, ipf) G S 2 } and a bounded set of controls (bounded for 
the norm C°°) {e(ipi,ipf) G ^t^^/AW; G 5 2 } such that for any two states ^ and ^ in S the 
solution of ([2]) with the control e(ipi,ipf) and starting from the initial state ^(0) = ipi arrives in ipf at 
time T (ipi, ipf). Hence, T — r 2 can be chosen bounded for all £. Therefore ^—U(0, r x ) and q^-U(t 2 , T) 
are bounded. Thus, we have: 

P l} ie) = 2M((f\U(T,T 2 ) — U(T 2 ,T 1 )U(T 1 ,0)\i)(i\U\T 1 ,0)U^T 2 ,T 1 )UHT,T 2 )\f)) + O(l) (16) 

°t l lk °t L lk 

We now utilize U(ti,0) \i) = \ipi) and W(T,t 2 ) \ f) = \ip 2 ) where \ipi) and |t/> 2 ) are defined above, and 
replace -A-U(ti, t 2 ) by its expression in ( fT5] ) to find: 



a /4 



g 
3/4 



2£ 



+ 0(1). 
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This expression holds for the control defined as: 

f ei(r), if re [0,n] 

<r)= I ^cos^r-n)), ifre]r 1; r 2 [ (17) 
I e 2 (r), ifre[r 2 ,T] 

z'z'J Part II: We now need to prove that -5-7— Pj/(e) = 0(1) for {m, n} 7^ {/, A;}, where e is the control 
found above in ([17]). As in equation < \l 1[ ), we have: 



^_/7(r 2 ,r 1 ) = z^(r 2 ,r 1 ) / cos(u4(r - n))C/t(r, n)<T™U(T, n)dr , (18) 



and again the result of lemma [2] is employed. Equation ( fT8| ) calls for 



cos(u4(r - r 1 ))e lH '^ T - T l) (7 mn e - l ^(r-r = 

i cos((u4 - ^J(r - - i sin((a/ ft - <J(r - n))a™ (19) 

+1 cos(K fc + <J(r - n))ar + 5 Bin((^ fc + ^J(r - n ))a^ . 
Considering that ff has non-degenerate transitions (see definition in appendix [A} implies that cu' lk — 



u' mn 7^ and u' lk + u' mn 7^ 0. As the expression in ( fT9[ ) oscillates at frequencies independent of £, it 



therefore contributes to 0(0 in ( [T8] ). Hence, for r 2 — ri = ^ we can directly conclude that: 





-Pif(e) = 0(1) . 



C. Proof of lemma [2] 

Proof: This proof relies on three consecutive changes of frame that aim to cancel the oscillating 
terms of order 1 and £ in the dynamics. We then derive a specific form of the averaging Theorem (see 
theorem 4.3.6 in [ 25 ] for a general form of the averaging theorem). 

For the sake of clarity and with no loss of generality, we take T\ = and note U(r) = U(t,t\). 
Equation ( flO] ) may be written in the interaction frame Ui(t) = e lH ° T U(r), 



where: 



^H t (t) = I /4n^ +a)U)T \m) <n| (20) 

(m,n)^(fc,Z) 

+ ^ E ^n^ + ^ Jr \m) (n\ 
and the average of Hj is zero. The average of a time dependent operator O(r) is defined as follows (see 



defintion 4.2.4 in (25)): 

0= lim ^ / 0(r)dr (21) 
0— >+oo y 







We now take Uj(t) = (I — i£Hi(t))Uj(t). Since -§pHj is almost periodic (cf. appendix |a| for the 
definition), then Hj is also almost periodic and hence bounded for all r. Hence, there exists £ > 0, 
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V£ < £ . and we have that / - z£#i(t) has an inverse and (I - i^Hfa))- 1 = I + i^H^t) + 0(£ 2 ). We 
find: 

= • U^a? - < 2 (^[^(r), a?] + Hfa^Hfa)) + O^)) Ufa) 



^[Hfa)^ 1 *] + H i {t)£H i {t) = i(K + ^(r)) where K = -%H 1 ^H 1 does not depe nd on j, and 
A'(t) is almost periodic of average zero and bounded for all r. It is important to note that \§^Hf = = 
Hj£Hj + {^HrjlTi = i(K - tf). Hence K = KUs Hermitian. 

We now take U"{r) = (I — i£ 2 K(t))U'j(t). Since K{r) is bounded for all r, then for a sufficiently 
small f, I-ii 2 k{r) has an inverse and (7-^ 2 if (t)) _1 = I + i£ 2 K(r) + 0(£ 4 ). [/}' satisfies the following 
equation: 

^ ^'(r) = i (^i k + i 2 K + 0(e)) U'{{t) , (22) 



<9r 

and we define U av to be the solution to the averaged dynamics: 



8 U av (r) = i (^a l x k + eK) l ir (r) (23) 



Or ^ ' ^ 2 

f/«,(0) 



We can directly solve p3[ ): 



Subtracting ([22]) from ([23]), we find, using Gronwall's lemma, that for all r < i one has U"(r) = 
U av (r) + Also note m at to go from [// to we have used two consecutive changes of variables 

which are close to the identity, hence: Vr U"(r) = Ui(t) + 0(£). Finally, since e~ tH ° T is an isometry, 
we have: , 

U(t) = e-^T^^+eVr + for all r < 

4 



V. Conclusion 

We have proved that the identification of the real dipole moment matrix of a controllable finite 
dimensional quantum system with non-degenerate transitions using as measurements only one population 
at a final time T is a well posed problem. That is, we can always find a control which discriminates 
between two close but different dipole moment matrices leading to two different measurements and this 
distinction is of first order. What at first may appear to be very restrictive data can actually be turned into 
a rich source of information by adequately controlling the dynamics of the system. 
It would be interesting to extend this result to the more general case of a Hermitian observable O besides 
the population projection operator. Also, the extension to the case where the dipole moment \i is an arbitrary 
Hermitian operator (not necessarily real with a nul diagonal) would be another important step. Finally, 
further directions include establishing under what conditions the identification is globally well posed (i.e, 
not just locally) and exploring the impact of laboratory noise. We are now investigating the experimental 



implementation of a Hamiltonian identification procedure along the lines of the proposal given in [26|, 
where the aim is to find, in the laboratory, those good sets of controls {ei, ..,e m } complying with all 
experimental constraints which make the inversion procedure well posed. 
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Appendix A 
Definitions 



In this appendix, we give some useful definitions: 

Definition 1. We say that system ([2]) is controllable [27] if and only if for all , \ift 2 ) £ <S there exists 
a time t and a control e G S t such that for \i[>(0)) = \4>i) © leads to \4>(t)) = \ip 2 ). 

Definition 2. We define a function J : fx G M. — >■ J{fi) G R. Let a > 0. Suppose J is C 2 . Let V 2 </(//.) 



be the Hessian of J at ft. We say that J is locally a-convex j\28)l around fi if the smallest eigenvalue of 
V 2 J(^i) is larger than a. 

Definition 3. Let Ho be a real N x N diagonal matrix with entries Ei,..,En. We say that Hq has 



non- degenerate transitions [29] ifV(l, k) ^ (m, n) I ^ k and m ^ n, we have \Ei — E k \ ^ \E m — E n \. 

Definition 4. Take system Let us denote M. as the space to which \i belongs. We say that \i is locally 
observable in M. if there exists r > such that for all fi G M. with < ||/2 — lm\\ < r there exists T > 
and e G St such that Pif(e,fi) ^ Pj/(e,/i). 



For more details on almost periodic functions, one can refer to [ 30 1 . Here, we give a very simplistic 
definition which is suitable for our needs. 

Definition 5. We say that an operator A(t) is almost periodic if there exists M G N and M time 
independent operators A , .., A M and M frequencies Wi,..,um such that for all t: A(t) = J2k=i e iu)kt A k . 

Appendix B 
Calculations 

A. Calculation 1 



We desire to compute 5U(t 2 ,ti) [31 1 with only ji varying, and all other elements in the dynamics are 



fixed. Taking the differential of equation ([5]), we get 

^-6U{r 2 ,Tt) = --^—{(H - e(r 2 )p)5U(r 2 ,r 1 ) - e(r 2 )6fiU(r 2 , n)) . (24) 

OT 2 LH 



We seek a solution of ( |24| ) in the form SU(t 2 ,ti) = U(r 2 , Ti)/3(r 2 ) where (3 is a matrix to be found 



satisfying /3(ti) = 0. Substituting this in (|24j), we find: 



1 



T-2 



P{t 2 ) = i— / e(r)U\r, n)5LiU(r, n)dr 



Hence, we have: 

r>T2 

— JKt^.tA I 

,\Ho 

and in particular: 



SU(t 2 , n) = 17^U(t 2 , n) / e(r)t7 t ( r, Ti)5liU{t, n)dr 



_d 
dfi 



rU(r 2 , n) = ijj^U{r 2 , n) / " e(r)U\r, n)a l x k U(T, n)dr (25) 

Ik \\ n 0\\ Jt! 
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B. Calculation 2 

In this section, we compute e~ t ^~ (7x + ^ K ^ T ^ T ^>a l ^e ll ^~ a ' x + ^ K )( T ~ T v. Since we assumed ji' lk ^ 0, 



A ^ a lk +i 2 K){T _ Ti) ^ k ^ { ^ a lk +eK){T _ Ti) 



e M 'fc (Tj, e "j* 



cr' fc + £^ is Hermitian, and is therefore diagonalizable. Hence, there exists a unitary matrix Pc and a 
real diagonal matrix A € such that <rj* + = P 5 x A 5 x Pj. The function £ G [0,f ] -> ^* + ^ is 



analytic, therefore the eigenvectors of cr^f + can be continued analytically as a function of £ (see 



A<!fc 



Theorem 6.1 in chapter II, §6 section 1 in [32|). Hence, we can write P^ = P + 0(£) where P is suc h 



that P x (jj x P T = o^ fe is real and diagonal. a l z k = \l) (l\ — \k) {k\. We have, Vr 



e (jj. e 



Poe _i^ e A € (T-n)pt <7 »p oe i>CAe(T-r l )pt + 
p oe -i^€A € (r-n) (T lfc e i4 fc €A e (r-r l )pt + (£) 

Po^pJ + oce) 



